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SUMMARY 


The  response  of  KC-135  and  747  type  aircraft  to  a  nuclear  burst  was  calculated 
using  the  VIBRA-4  code.  Runs  were  made  with  and  without  the  diffraction 
phase  loading  to  determine  its  effect  on  the  aircraft  response. 

Bending  moments  were  computed  for  stations  on  the  fuselage,  wing,  horizontal 
and  vertical  stabilizers,  for  bursts  from  below,  above,  and  to  the  side.  The 
ratios  of  peak  bending  moments,  with  and  without  the  diffraction  phase  were 
calculated.  The  differences  in  the  calculated  ratios  were  small  for  the 
KC-135  model,  but  large  differences  were  found  for  the  747  models.  It  is 
clear  that  the  diffraction  phase  loads  can  exert  a  large  influence  on  the 
nuclear  gust  imposed  loads,  but  it  is  not  possible  to  predict  this  effect 
with  any  confidence  in  advance  of  detailed  calculations. 

Thus  any  computer  codes  for  nuclear  gust  response  should  include  the  diffraction 
phase  loads.  Further  work  should  be  made  to  verify  loading  trends  identified 
in  this  study. 


CONVERSION  FACTORS  FOR  U.S.  CUSTOMARY 
TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


To  Convert  From 


To 


Multiply  By 


foot 

foot/ second  (fps) 

Inch 
kllotons 
kip  (1000  lbf) 
klp/lnch*  (ksl) 
knots 

nautical  miles 

pound-force 
(lbf  avoirdupois) 

pound- force  Inch 

pound-force/ Inch  (psi) 

pound-mass 
(Ibm  avoirdupois) 

2 

pound-mass -foot 
(moment  of  Inertia) 

pound/square  foot  (psf) 
slug 

square  foot  9 
(sq.  ft.,  ft*) 


meter  (m) 

3.048  000  X  E  -1 

meter  per  second 

3.048  000  X  E  -1 

(m/s) 

meter  (m) 

2.540  000  X  E  -2 

tera joules 

4.183 

newton  (N) 

4.448  222  X  E  +3 

kilo  pascal  (kP$) 

6.894  757  X  E  +3 

meter  per  second 

5.144  444  X  E  -1 

(m/s) 

meters  (m) 

1.852  000  X  E  +3 
• 

newton  (N) 

4.448  222 

newton-meter  (N  m) 

1.129  848  X  E  -1 

kilo  pascal  (kPa) 

6.894  757 

kilogram  (kg) 

4.535  924  X  E  -1 

2 

kilogram-meter 
(kg  mz) 

4.214  011  X  E  -2 

2 

kilogram  per  meter 
(kg/mz) 

4.882  428 

kilogram  (kg) 

1.459  390  X  E  +1 

meter*  (m*) 

9.290  304  X  E  -2 

A  more  complete  listing  of  conversions  may  be  found  In  "Metric 
Practice  Guide  E  380-74,"  American  Society  for  Testing  and 
Materials. 
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1.0  INTROnUfTIDN 

In  order  for  our  strategic  forces  to  remain  a  credible  deterrent,  i  +  's 
necessary  tha*  they  have  a  demonstrated  capability  to  perform  their  m;ss>or 
This  demonstration  has  many  aspects  and  one  of  these  is  to  show  *\v  an^'vc  s  fhai 
sufficient  weapon  systems  will  survive  an  attack  and  he  able  to  completp  thp-r 
mission  of  inflicting  unacceptable  damage  on  the  attacker.  Pol ’rv  makers  musr 
have  high  confidence  in  such  analyses  for  fhe  analyses  will  revea1  the  strengths 
and  weaknesses  in  our  various  weapon  systems.  Then  the  policy  makers  will  guide 
the  development  of  improvements  or  change  in  our  defense  system.  Te  tha+  end 
the  following  study  was  conducted  to  evaluate  a  part  in  ’ dr  portion  o*  fhr 
analytic  tools  available.  This  study  was  concerned  wifh  *np  impertdnre  o* 
modeling  the  aerodynamic  forces  on  an  aircraft  in  the  time  period  irrmed'ately 
after  shock  arrival. 

The  aerodynamic  forces  acting  on  an  airfoil  are  the  result  of  pressure  on  the 
airfoil  surfaces.  The  difference  in  pressure  between  upper  and  lower  surfa- es 
of  a  wing  for  example  results  in  1’ft.  <Vpssure  depend'  on  the  a'rf’ow  uattern 
around  the  airfoil.  When  the  airflow  pattern  is  abruptly  changed  t  tip  re  is  a 
time  of  transition  until  a  new  airflow  pattern  and  new  pressure  patterns  arp 
est.abl  ished. 

A  shock  wave  passing  over  an  airfoil  1  or  any  structure'  produces  unbalanced 
overpressure  forces.  initially  oril v  the  surface  first  encountered  experiences 
overpressure,  then  more  and  more  of  the  surface  experiences  nvprp>' “ssurp  as  the 
shock  wave  engulfs  the  structure.  At  the  initial  encounter  the  shoe*  wave  ’ c 
reflected  and  the  front  surface  experiences  the  reflpctpd  shock  over pressurp. 
When  the  shock  wave  completely  engulfs  thp  structure  a  rarefar+’on  wave 
propagates  back  to  the  front  face  the  structure  and  reduces  the  rp*ie,  epo 
shock  pressure  to  its  unrr fleeted  value.  This  period  is  often  referrpd  *<■  a c 
the  diffraction  phase.  Aerodynamic  loads  during  fhe  differ*  on  phase  ire 
uncertain  for  two  reasons.  First,  the  proper  reflexion  factor*;  ♦  o  i se  for  <-br 
complex  aerodynamic  shape  of  airfoils  are  not  accurately  known.  Measurements  of 
reflection  factors  have  been  made  for  simple  shapes,  such  a*  ’nrlinpd  rl  anps  and 
cylinders,  but  not  for  complex  shapes.  Second,  the  time  for  aerodynamic  lift  to 
change  after  an  abrupt  flow  change  is  not  accurately  known.  It  was  not  known 
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whether  diffraction  phase  loads  were  dominant  in  causing  damage  to  an 
aircraft  in  a  nuclear  encounter.  If  they  dominated  over  later  aerodynamic 
flow  loads,  then  better  resolution  of  diffraction  phase  uncertainties  would 
be  necessary.  If  they  were  not  dominant,  then  present  understanding  would  be 
adequate.  An  investigation  of  the  importance  of  the  diffraction  phase  loads 
was  made  by  changing  the  diffraction  phase  loads  in  a  nuclear  gust  response 
program  and  determining  the  resulting  effects  on  structural  response. 

The  structural  response  was  calculated  for  two  different  aircraft  using  the 
VI8RA-4  code  with  and  without  the  diffraction  phase  loads.  The  results 
show  that  although  the  diffraction  phase  loads  are  not  large  enough  to  cause 
major  structural  damage  it  is  important  to  include  them  with  the  gust 
loads  to  determine  the  structural  response.  For  some  vehicles  at  critical 
locations  the  peak  bending  moment  calculated  by  VIBRA-4  with  the  diffraction 
phase  loads  included  were  two  to  three  times  as  large  as  those  calculated 
without  them. 
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2.0 


COMPARISON  OF  DIFFRACTION  PHASF  MF THORS 


Several  methods  have  been  developed  and  documented  to  handle  loads  on  air 
vehicles  due  to  nuclear  blast.  A  few  of  the  methods  are  described  briefly 
below.  Many  more  methods  and/or  computer  programs  exist  with’"  ^he  many  orivate 
companies  that  are  concerned  with  the  effects  on  air  vehicles  of  nuclear  blast. 

?.l  DNA  2D48H-1 

A  general  description  of  structural  analysis  for  qtist  response  and  for 
diffraction  phase  response  is  described  in  Section  III-7.3  of  the  Handbook  for 
Analysis  of  Nuclear  Weapon  Effects  on  Aircraft  (Reference  7).  in  this  reference 
"overpressure  effects"  describes  the  same  loading  which  is  called  diffraction 
phase  loading  in  this  document. 

In  general  this  reference  separates  the  effects  of  gust  response  and  the  effects 
of  the  diffraction  phase  loads.  Because  of  the  duration  of  each  type  0f  loading 
and  the  structural  damage  caused  by  each  type  of  loading,  *he  analyst  can 
consider  them  separately.  For  the  hand  calculations  described  and  the 
accuracies  which  these  methods  yield  the  above  assumption  may  be  a  valid 
assumption. 

2.2  NOVA-? 

The  NOVA-2  digital  computer  program  (Reference  4)  is  an  extension  of  the  nova 
program,  and  is  one  of  the  Method  ’  calculations  described  m  dna  2048H-1.  The 
NOVA  and  NOVA-2  programs  are  digital  programs  to  calculate  the  linear  or 
non-linear  response  of  localized  aircraft  structures  to  overpressure  or 
diffraction  phase  loads.  For  these  programs  the  gust  response  time  is  long 
compared  to  the  structural  response  of  interest  and  therefore  is  not  included  m 
the  program. 

2.3  VIBRA-4 

The  VIBRA-4  digital  computer  program  (Reference  1)  is  one  of  the  Method  • 
calculations  described  in  DNA  2048H-1.  This  program  is  primar’ly  a  gust 


response  program.  However  the  aerodynamic  forces  include  an  abrupt  change  due 
to  shock  arrival  and  passage  over  a  surface.  The  rest  of  the  computation  time 
is  for  gust  response.  This  is  the  program  selected  for  this  study.  The  VIRRA-a 
program  used  for  this  study  is  more  fully  described  in  Section  1.0. 

2.4  VIBRA-6 

The  VIBRA-6  digital  computer  program  (Reference  S)  is  an  extension  of  the 
VIRRA-4  methodology.  VIBRA-6  is  a  new  program  which  uses  doublet  tattice 
aerodynamic  theory.  The  program  makes  the  calculations  in  the  frequency  domain 
and  uses  a  Fourier  transform  to  generate  the  time  histories  of  the  gust 
response.  In  Section  1 1  -4  of  Volume  II  (Reference  b  >  n*  the  VIRRA-6 
dor umentat ion  it.  states  that  VTRRA-6  currently  does  not  contain  +  he  equations  to 
include  the  effects  of  the  diffraction  phase  loading.  the  resists  shown  in 
Section  4.0  of  this  study  show  that  the  equations  in  Appendix  A  of  Reference  6 
should  be  added  to  the  VIBRA-6  program. 
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GUST  ANALYSIS  METHODOLOGY 


3.0 

Thp  gust  environment  associated  with  a  nuclear  burst  is  manifested  by  a  sudden 
change  in  local  density  and  relative  wind  velocity  as  fhe  shock  wave  sweeps 
across  the  aircraft.  These  density  and  velocity  changes  <-psu1  t  m  two  d'st’nc' 
types  of  loading  on  the  aircraft  lifting  surfaces  and  fusel aqe. 

The  first  is  called  the  diffraction  loading  and  ’asts  approx imate1y  tne  time 
’-equired  for  the  shock  wave  to  travel  one  chord  of  the  aerodynamic  surface  beino 
considered.  Durinq  this  time,  the  flow  is  +wo-dimensional  and  the  normal  force 
coefficient  is  given  by  rn  =  4/M,  where  M  is  the  free  stream  Mach  number  based 
on  conditions  behind  the  shock  fron*  . 

The  second  part  of  the  loading  v  obtained  by  considering  fhe  three-dimensional 
lift  buildup  due  t0  circulation  >ffec*s.  During  rapid  ’oading  provided  by  the 
ma'enal  velocity,  the  lift  rirculation  follows  a  Do’tdup  given  by  a 
re’ationshiD  similar  to  the  Wagner  function.  This  buildup  m  lift  is  actually  a 
family  of  curves  depending  upon  fhe  annle  f  attack  of  the  local  lifting  sermon 
and  the  chord.  Since  +hp  material  velocity  effects  all  the  aerodynamic  surfaces 
almost  simultaneously,  analysis  of  an  a’ '"craft,  in  this  environment  must  consider 
the  complete  aircraft.  Another  reason  why  the  analysis  mus+  consider  the 
romplpte  craft  is  that  the  gust  may  enqulf  one  part,  fe.g.,  right  wing  tipi  and 
affect  the  position  and  orientation  of  other  points  (°.g.,  left  wing  tipi  beforp 
the  shock  arrives  at  that  point.  Similarly,  loading  subsequent  to  the  arrival 
of  the  shock  wave  depends  upon  the  elasticity  of  the  aircraft  arKj  interaction  of 
each  major  component  upon  the  other. 

Therefore,  to  assess  the  vul nerabi1  ity  of  the  aircraft  in  a  gust  environment,  a 
dynamic  model  must  be  developed  which  is  actually  a  mathematical  description 
representing  the  eiastic  and  inertia  properties  of  the  aircraft  once  transient 
aerodynamic  loads  are  applied  to  the  lifting  surfaces.  In  addition,  the 
analysis  must  be  achieved  using  computerized  techniques  if  Hmelv  analyses  are 
to  be  conducted.  The  computeri zed  technique  used  in  this  analysis  is  described 
in  Reference  1.  The  method  is  given  thp  name  VIBRA-4,  and  althouqh  it  is 
described  in  detai1  in  the  cited  references,  a  brief  description  of  the 
methodology  will  be  given  to  provide  an  understanding  and  eva’uation  of  changes 
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to  the  technique  which  were  made  during  the  coursp  of  this  analysis.  The 
equations  of  motion  are  described  first,  followed  by  a  discussion  of  i-he 
aerodynamic  loading  methodology  and  finally,  a  discussion  of  failure  criterion. 

3.1  V 1  BRA EQUATIONS  OF  MOTION 

The  equations  of  motion  used  in  VI8RA-4  are  derived  from  Laqranqes  equation, 
which  is  given  by 


_d  /  3T  \  3J  JV 

dt  \  3q  /  3q  3q 


(1) 


is  the  kinetic  energy  of  the  system 
is  the  potential  energy  of  the  system 
is  any  generalized  coordinate  of  the  system 
is  the  generalized  force  associated  with  the  coordinate  q 

The  kinetic  energy  term  will  be  described  first,  followed  by  the  potential 
energy  term  and  the  generalized  force. 

3.1.]  Kinetic  Energy 

The  kinetic  energy  term  in  Fquation  1  is  nf  the  form 


where 

T 

V 

q 


Tl  /p  \ 

T  =  1/2  nuv.2 

i=l 

where  ^ 

mi  is  the  i  mass  of  the  system  of  discrete  masses  describing 

the  mass  distribution  of  the  system 

v.j  is  the  linear  velocity  of  the  i**1  mass  point  of  the  system 

n  is  the  number  of  masses  describing  the  complete  system 

T  is  the  kinetic  energy  of  the  system  with  respect  to  a 

reference  point,  specifically  the  center  of  gravity 
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where 


V  y0  and  }Q  are  thp  components  o*  accel  era*  i  on  o*  thp 

rpfere nee  point  in  *he  FAAS 

i-  ,  1-  and  1-  are  thp  components  of  acceleration  of  mass  point 
i  yi  i  i  relative  to  the  reference  point  in  the  FAAC 

The  values  of  £;■  ,  i~  and  IL'j  are  obtained  in  the  following  manner.  The 

xi  y  i  zi  , 

accel erat i ons  of  mass  points  will  be  calculated  in  the  moving  aircraft  system, 

therefore,  a  coordinate  transformation  from  the  MAAS  system  to  the  FAAS  is 

required.  Using  the  coordinates  of  ‘he  mass  point  in  the  FAAS  and  the  MAAS  to 

be  ,  £p  ,  £?  and  %  ,  £  ,  and  £  ,  respectively;  and  applvinq  an  Filler 

xi  yi  zi  xi  yi  zi 

transformaf ion  the  following  is  obtained 

l-  =  £x  (cos  0  cos  £)  +  £  (sin  8  sin  0  cos  j>  +  cos  5  sin  -Ji) 

i  yi 

+  i  (-  cos  C  sin  $  cos  $  +  sin  B  sin  5) 

1'  ffii 

£.-  =  £  (-  cos  §  sin  ijj)  +  £  (cos  6  cos  ij'  -  sin  B  sin  9  sin  tb ) 

vi  xi  yi 

+  £z  (sin  B  cos  $  +  cos  8  sin  9  sin  >£) 

l-  =  t  sin  O  -  £  (sin  B  cos  0)  +  £  (cos  B  cos  0) 

Zi  xi  yi  zi 

where 

e,  B,  are  the  rigid  body  rotation  angles  o*  the  reference  point 
in  the  FAAS. 

Equation  6  is  then  differentiated  twice  with  respect  to  f’tne  *o  obtain  the 
relationship  between  the  relative  accel erat ions  in  +he  FAAs  ,nd  the  MAAS.  The 
next  step  is  to  obtain  the  rotational  accel erations  "n  +  hp  caas  in  terms  of 
rotations  in  the  MAAS.  These  rel at ionshi ps  become 

8-6  -*9  (7) 

9  =  0  +  ipB 

ip  =  ij;  -  96 

Making  these  substitutions  into  the  acceleration  equations  resulting  from  the 
differentiation  of  Equation  6  and  letting  the  FAAS  and  the  MAAS  be  coincident  at 
any  instant  of  time,  a  substantial  simplification  in  the  equations  (i.e.,  all 
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7  ‘ 


WK  / 


cosine  terms  become  1.0  and  all  sine  terms  become  0.0)  can  be  made. 
The  final  equation  becomes 


£  =  i  +  £  -  l  e  +  21.,  hi  -  2i  e  -  £  (e2  +  <P)  +  b(  £  «  +  £  i) 


xi  xi  yi  zi  yi 


£  =  1  -  £  +  £  6  -  2£y  +  2i  B  -  l  (62  +  J2)  +  e(  £x  B  +  £z  i) 


yi  Yi  X.  z.  X.  ^  y1 


£  =  ii  +  £  0  -  £  8  +  2£y  e  -  2i  6  -  £_  (B2  +  62)  +  <!>(£  B  +  £  0) 

zi  zi  xi  yi  Xi  Yi  Zi  Xi  yi 


One  final  step  is  taken  as  part  of  the  norma1  mode  approach.  Let  the  coordinate 
of  a  mass  point  at  any  instant  of  time  be  the  original  coordinate  plus  its  modal 
displacement  as  shown  in  Equation  d 


£  =  £  +  V  <J>  ^  q 

xi  xi  kl  Xi  r 


£  =1  +  V  <t>  ^  q 

yi  yi„  r?i  r 


■  =  £  ^  <(,  (»*) 

*z.  Xzi  +  £  %  qr 

1  o  ’ 


where 


£x  ,  £y  .£z  are  the  coordinates  of  mass  i  when  the  aircraft  <s 
i  i  undeformed 


ooo 


is  the  nunber  of  elastic  modes 


i  is  the  deflection  of  the  rx  mode 

nr 

♦  ^  are  the  modal  deflections  in  the  rth  mode  of 

xi  yi  zi  mass  point  i 
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The  first  and  second  time  derivatives  of  these  coordinates  then  become 


and 


"E 

£ 

r=l 


"E 

£ 

r=l 


nE 

£ 

r=  1 


"E 

£ 

r=l 


"E 

£ 

r=l 


"E 

£ 

r=l 


M  O  i 


(in 


A  similar  procedure  is  used  to  relate  the  rigid  body  translational  accelerations 

in  the  FAAS  and  MAAS.  The  resulting  equations  are 

x  =  x  +  iiiy  '9z 
o  o  o  o 


y  =  y  +  Bz  ~  \|)  x 

o  o  o  o 

•  •  •  •  •  • 

z  ~  z  -  By  +  0  x 

o  o  o  o 


(!?' 


Equation  12  is  found  by  differentiating  the  relationship  between  velocities  in 
the  FAAS  and  the  MAAS,  which  are  similar  in  form  to  those  of  Equation  6,  and 
applying  the  simplification  of  assuming  that  the  FAAS  and  MAAS  are  coincident  at 
any  instant  of  time.  One  other  set  of  equations  is  needed  to  support  Equation 
12.  This  is  the  set  of  equations  relating  the  rigid  body  rotational  velocities 
in  the  fixed  and  moving  systems.  These  equations  become 

.  t  - 

B  =  8  +  4>  sin  0 

9  =  e  cos  8  -  I  cos  9  sin  8  H’ 

#  f,  •  —  m 

=  0  sin  B  +  $  cos  0  cos  8 


3.1.?  Potential  Energy 


The  potenti a1  energy  of  the  system  is  composed  of  the  potential  energy 
contributed  by  the  elastic  vibratory  modes  of  the  structure  and  that  resulting 
from  gravitational  forces.  The  contribution  from  the  elastic  vibratory  modes  of 
the  structure  may  be  written  as 


wherp 


M 


r 


Vv  •  1/2  £ 

r 


M  a)  * 
r  r 


f  14) 


is  the  general i?ed  mass  associated  with  thp  rth  norma1  mode 
of  the  system 


% 


k(r)  1 

2 

+ 

[♦  (r)l 

2 

+ 

- 1 

e- 

-$ 

2i 

L  1 

j 

J 

vibrational  freque^^y  0f  the  normal  mode 

is  the  rf^  normal  coordinate 

represents  the  r  normal  mode  of  the  system 


Equation  15  gives  the  contribution  from  gravitational  forces 


V 


9 


sin 


'151 


sin  Be 


cos  0£ 


/-  V 

—  • 

1  z  +£,,  \ 

cos  Br  cos 

\  0  zi  / 

E  E 

where 


®E»  arp  f:he  £u1pr  angles  relating  the  coincident  axes  *o  the 


Earth  Fixed  Axis 

m.  represents  the  i  mass  of  the  system 

The  total  potential  energy  then  becomes 

V  =  V  +  V 

v  g 


(161 
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3.1.3  Generalized  Forres 


The  generalized  force  associated  with  t.hp  coordinate  q  may  he  represented  by 


q  -|SL 

\  3qr 


where 


represents  the  virtual  work  of  the  system 
is  the  rt^1  generalized  coordinate 


The  partial  derivative  of  the  virtual  work  is  given  hy 


3W_  _  F _ I 

'  x  3q, 


i  +  F  yi  +  F  "i 


where 


F  F  F  are  the  external  forces  /aerodynamic,  etc.'  acting  on  the 
x  ’  y  ’  z 

system 

3.1.4  Equation  of  Motion 

Using  the  components  presented  in  the  previous  paragraphs  and  conducting  the 
operations  given  by  Equation  ,  "he  equation  of  motion  used  in  VIPRA-4  is 
obtained. 

£  j  \<r)  i +  ♦», (r>  51  +  %  w  51 )  *  w2  ^ 


+  m^g  4>x(r^  sin  eg  -  <i>y^  sin  BE  cos  §E  +  <J>2(r)  cos  BE  cos  0£ 


V  I  (’ 

i*l  \  xi  xi 


yi  *i  zi  zi 
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where  the  components  of  Equation  19  are  previously  denned.  rouation  ’9  and  ”s 
supporting  equations  may  be  written  as 

rMrs'!  (qs)  =  (RHSrl  '?n\ 

where 

RHSr  (righthand  sidel  is  a  function  of  velocities  and 

displacements  and  may  be  found  by  the  previously  defined 
equation 

q‘s  is  the  modeal  acceleration  components 

M  is  the  general i zed  mass  matrix 

The  solution  procedure  is  accomplished  by  solving  for  the  acceleration  of  the 
generalized  coordinates  and  integrating  these  components  to  determine  velocities 
and  displacements.  The  integration  is  conducted  numerically  by  a  fi^th-order 
Adams  integration  formula. 

?.?  AFROPYNAMTC  LnAniNP 

The  forces  Fx  ,  Fy  and  Fz  given  in  equation  19  are  comprised  of  *wo  aistinrf 
sources  of  aerodynamic  loadings.  The  first  is  the  lg  flight  loads  which  are 
normally  acting  on  the  aircraft  and  the  second  is  the  incremental  loadinq 

produced  by  changes  in  material  velocity  and  density  behind  the  shock  wave.  The 
VIBRA-4  program  was  modified  in  this  area  of  steady  loadinq.  These  changes 

(made  by  AFWL)  will  be  described  below.  This  version  of  the  VIRPA-4  program  was 
restricted  to  subsonic  aerodynamics  only.  The  previous  versions  of  the  program 
did  not  consider  a  variation  in  the  ratio  of  lift  coefficient  to  angle  of  at+ack 
across  the  span  of  components  and  it  assumed  that  the  center  e*  pressure  was 
located  at  the  50  percent  chord  for  steady  loads  as  well  as  the  unsteady  load 
where  50  percent  is  appropriate.  Modifications  were  made  to  set  the  center  of 

pressure  at  the  40  percent  chord  which  was  a  comprise  between  the  lq  flight 

condition  and  the  unsteady  load  condition.  A  subsonic  load  subroutine 
(Weissinger  method)  documented  in  Reference  9  was  added  to  the  VTRRA-4  program 
to  provide  the  steady  airload  distribution  calculations. 

To  determine  the  lg  flight  loads,  the  program  adjusts  the  wing  and  horizontal 
stabilizer  angle  of  attack  to  provide  a  balanced  load  condition.  During  this 
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♦  *  \  f  i 


♦pr<r  On  D0r’Od,  *  hp  Strtictu^P  ^  ro  de'Hec*  ’  *  **  'ticxIh' 

a' low  ’t  .  At  fhp  end  of  rhp  Steady  fl  iqh*  calculation,  fhf  fnr<  o<-  ar’’nQ  or  "ip 
aircraft  are  completely  balanced  and  thp  loads  at  each  o+  'hr  dr " '  t'pr'  ,t'  'wt 

arp  ral ousted. 

The  sprond  component  of  aerodynamic  loading  results  from  thp  ma+er'al  velocity 
hpt’,nd  the  shock  wave.  the  mode1  for  calculating  *hp  overpressure  and  thp 
associated  material  velocity  components  for  a  given  vield  and  cange  is  discussed 
in  detail  in  Section  ‘>.7.  Thp  material  velocity  time  history  induces  a 
time- dependent  angle  of  attack  which  causes  time-dependent  load  to  ex  is*  'n  *he 
aircraft  striK  ture.  The  mstan* aneo.is  angle  of  attack  ,c  laicuiated  by 
assigning  thrpp  of  the  discrete  masses,  pither  within  or  near  an  aerodynamic 
strip,  to  pach  of  thp  aerodynamic  panel  <  so  that  a  plane  passed  through  the 
masses  represents  the  orientation  of  a  part-cular  strip.  The  unit  norma1  vector 
for  this  plane  is  then  found,  and  *he  anqle  of  attack  determined  by  forming 
dot  product  of  this  vector  with  *hp  velocity  vector  of  the  surfacp  relative  to 
the  air.  The  angle  is  measured  ’n  a  plane  containing  botn  the  velocity  and  the 
norma1  vector,  Fquation  ?\  shows  this  cel  at lonship. 

N  •  VA  =  |N||VA|cos  (90°  +  ,)  =  |N|  |VA|sin  n  <2\  ' 

where 

N  is  the  normal  vector 

VA  is  the  velocity  vector  of  thp  aerodynamic  strip  relative 
to  the  air 

a  is  the  angip  of  attack 

flnce  the  angle  of  attack  or  change  in  angle  of  attack  <  o  -  *0!  s  found,  the 
normal  force  coefficient  can  then  be  determined.  Subsequent v,  thp  li*t  force 
is  determined  from  Equation  77 . 


lA 


1  /?  o  V, 
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I 


where 


Lfl  is  the  ♦’otal  lift  acting  on  a  lifting  segment 

o  is  the  density  of  the  air  flowing  over  the  surface 

is  the  velocity  of  the  aerodynamic  strip  relative  to  air 
is  the  area  of  the  particular  strip 


The  components  of  the  normal  force  are  then  found  by  mul  t .i pi  yi  nq  the  tota1  force 
by  the  components  of  the  unit  normal  vector.  The  components  are  given  by 


N 

x  "  La  TnT 

Ny 

F  =  |  _ L 

y  A  |N| 
Nz 

Fz  LA  ]NJ 


r>->  \ 


These  forces  are  then  distributed  among  the  individual  masses  and  the  oespons° 
loop  continues. 

The  normal  force  coefficients  used  ace  from  two  sources,  depending  on  whe+hpr 
the  trim  or  steady  load  calculation  is  being  performed,  or  the  unsteady  loading 
due  to  gust  is  being  performed.  As  discussed  previously,  the  steady  load 
calcu1  at.ions  are  achieved  by  using  the  Weissinger  method.  For  the  unsteady 
portion  of  the  loading,  experimental  normal  force  coefficient  da+a  make  up  the 
portion  of  the  time  when  circulation  lift  is  predominant,  i.e.,  after  the 
aerodynamic  panel  has  been  completely  engulfed  by  the  shock  wave.  Another  wav 
to  indicate  this  time,  and  the  one  used  in  VIRRA-4,  is  by  using  the  concept  of 
aerodynamic  time  S  =  Vt/C,  where  V  is  the  vector  sum  of  the  velocity  of  the 
lifting  surface  and  the  shock  front,  t  is  time,  and  C  is  the  structural  chord  of 
the  lifting  panel.  When  S  =  0,  the  shock  has  just  arrived  at  the  lifting 
surface  panel.  When  S  =  1,  the  panel  has  been  engulfed  by  t,bp  shock.  Figure  0 
shows  the  basic  two-dimensional  unsteady  normal  forces  coefficient  to  angle  of 
attack  ratio  as  a  function  of  aerodynamic  time.  These  data  werp  obtained  in 
shock  tube  experiments. 

In  reviewing  the  unsteady  aerodynamics  of  the  VIBRA-4  program,  it  was  observed 
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that  each  lifting  surface  and  lifting  surface  panel  was  given  '-he  sane  norma1 
force  coefficient  to  angle  of  attack  ratio;  i.e.,  the  value  of  f^/a  was  <onsran+ 
across  the  span  and  for  each  surface.  To  avoid  this  limitation,  results  of  the 
airload  distribution  formed  during  the  steady  load  calculations  were  saved  and 
used  to  scale  the  data  of  Figure  ?  to  be  more  applicable  to  ’he  idtinq 
surfaces.  The  net  effect  o*  this  scaling  process  is  to  reduce  the  tota'  loading 
slightly  and  to  redistribute  th  loading  across  the  surface  to  bettor  match  the 
actual  distr ibut ion . 


rusp1age  aerodynamic  loading  is  obtained  in  a  similar  manner  except  that  the 
force  is  nroportional  to  a  cross-flow  drag  force  given  by 


F 


N 


o 

0.3Sy 


1 


where 

F^  is  the  cross-flow  drag  force  in  the  direction  of  VN 

F 

o  is  the  ambient  air  density  at  the  fuse1 age  aerodynamic 

segment 

S  is  the  area  of  the  fuselage  aerodynamic  segment 

VN  is  the  flow  velocity  relative  to  the  aircraft  anH  normal  to 

T  the  fuselage  1ongitudinal  axis 


No  aerodynamic  loading  is  included  for  the  engine  nacelles.  Inclusion  of 
nacelle  aerodynamics  should  provide  damping  to  the  nacelle  motion,  hut  it  will 
also  provide  some  limited  nacelle  loading.  Jt  is  unknown  just  how  much  effect 
nacelle  aerodynamics  would  make  in  the  present  results.  more  ’n-depth  study 
should  look  into  this  particuiar  problem  and  include  a  morn  adequate  failure 
moment  calculation  using  combined  loads  instead  of  the  single-side  load 
presently  used. 

3.3  GUST  ENVIRONMENT 

A  new  AFWL  1  KT  nuclear  standard  (Reference  10)  based  on  data  from  the  AFW1 
SPUTTER  and  SAP  fluid  dynamics  programs  provides  the  free  air  overpressure, 
material  velocity,  and  overdensity  waveforms  for  the  blast  mode1.  This  nuclear 
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standard  replaces  the  modified  I RM  Problem-M  dafa  used  in  the  VIRRA-4  Mas* 
model.  Figure  4  indicates  that  for  overpressures  less  than  R  to  10  ps- ,  ’h-- 
standard  predicts  a  range  which  is  appreciably  shorter  than  that  obtained  us’rq 
the  !RM  Problcm-M  data. 

Replacing  the  free-air  problem-M  curve  with  the  AFWL  «  KT  sea  level  stannarc 
required  a  modification  to  two  scaled  reflection  factor  curves  given  in  Sect  on 
H  of  Peference  1.  These  factors  describe,  as  a  function  of  ranqe,  *he  blast 
r harar^eri st ics  in  thp  regions  of  regular  reflection  and  Mach  reflection.  Thp 
scaVd  reflection  factor  CRFRRR  is  based  on  height  of  burst  (HORl  curves  and 
applied  to  the  reflection  phase  of  a  near-ground  burst.  A  review  of  HOP 

overpressure  data  indicated  that  the  overpressure  along  the  ground  in  ’‘he  region 
of  regular  reflection  could  be  reproduced  effectively  by  considering  an 
imaginary  nuclear  explosion  with  a  specified  effective  yield  at  the  original 
burst  point.  The  effective  yield  can  be  obtained  by  modifying  the  actual  yield 
by  a  reflection  factor.  In  this  manner,  a  curve  of  the  cube  roof  of  this 
reflection  factor  for  the  regular  reflection  region  CRFRRR  was  determined  as  a 
function  of  the  slant  range  to  specific  points  on  thp  qround.  Likewise  the 

scaled  reflection  factor  CPFRMR  was  determined  for  thp  region  of  Mach 
reflection.  These  factors  are  used  in  the  blast  model  to  determine  the 

effective  yield  required  to  duplicate  the  overpressure  along  the  ground  in  the 
region  of  regiPar  and  Mach  reflection. 

1.4  FLIGHT  LOAD  ALLOWABLES 

Vulnerability  to  gust  effects  is  established  by  determining  the  burst  time/rangp 
at  which  either  a  limit  flight  load  (sure-safe)  or  ultimate  allowable  flight 
load  fsure-kill)  is  attained  at  some  location  in  the  aircraft  structure.  The 
ViBRA-4  code  calculates  all  forces  and  moments,  but  currently  provides  the 
capability  for  monitoring  only  up  and  down  bending  moments  on  the  wing, 

horizontal  stabilizer,  and  aft  fuselage  and  lateral  bending  moments  on  the 
vertical  tail,  aft  fuselage  and  nacelles.  Component  stations  having  thp  lowest 
design  margins  of  safety  were  selected  for  each  component. 

The  u’timate  allowables  were  calculated  by  modifying  the  ultimate  loads  by  thp 
margin  of  safety,  such  that  *he  margin  of  safety  is  redurpd  to  zero.  The 
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sure-safe  limit  allowables  in  this  table  are  two-thirds  of  the  ultimatp 
al lowables . 

Mission  completion  for  gust  response  was  considered  to  be  very  '■’ose  to  rde 
sure-kin  response;  i.e.,  if  the  load  is  just  slightly  less  than  ultimate  the 
structure  will  not  fracture  and  the  component  can  still  carry  loads  which  will 
allow  mission  completion.  This  result  is  substantiated  by  a  process  called 
"proof  test  logic"  which  effectively  guarantees  that  the  structure  withstands  a 
particular  load. 

Changes  made  to  the  VIRRA-4  code  to  more  accurately  trea+  romhined  vertical  and 
lateral  fuselage  loading  are  discussed  below.  The  program  changes  arrompl i shed 
to  calculate  internal  stresses  at  wing  and  fuselage  stations  were  disabled  *or 
the  current  study. 

The  allowables  for  wing,  horizontal  and  vertical  stabilizers  and  nacelles  are 
pure  bending  allowables.  Since  these  components,  except  for  the  nacelles,  act 
like  a  beam,  the  pure  bending  allowable  is  probably  adequate. 

The  fuselage  allowable,  however,  is  based  on  a  side-bending  moment,  an 
up- bending  moment  and  a  down- bend inq  moment.  The  original  VI BRA-4  methodology 
used  a  combined  loading  scheme  which  allowed  the  combined  loading  to  follow  the 
boundary  Rj  in  Figure  T.  This  allowed  both  the  ultimate  up  or  down  load  to 
exist  with  the  ultimate  side  loading.  For  pure  up  or  down  bending  or  side 
bending,  the  result  is  correct,  but  for  combined  loading  an  elliptical  boundary 
given  by  B2  in  Figure  ?  is  the  correct  solution.  Combined  loading  ’s  shown  for 
example  by  Pj  which  is  not  acceptable,  whereas,  Pp  s^ows  an  example  of  a  correct 
result.  The  ellipse  determined  from  a  combination  of  pure  side  bending  or 
vertical  bending  conditions  is  seen  to  provide  an  accurate  representation  of  the 
family  of  straight  line  segnents.  The  equations  used  in  the  program  for  the 
bending  moment  ratio  M/Mu-|t  are  given  by 


where 


Bu  is  the  allowable  up- bending  moment 

Bs  is  the  allowable  side-bending  moment 

Mz  is  the  Tabulated  up-bending  moment 

Mx  is  the  calculated  side-bending  moment 

R  is  the  ratio  of  calculated  moment  to  ultimate  moment 


This  equation  is  valid  for  up-bending  only.  If  the  bending  moment  is  down,  the 
correct  equation  is 


(26) 


where 

B^  is  the  allowable  down-bendinq  moment 

A  similar  technique  should  be  used  tor  the  engine  nacelles,  but  time  did  not 
permit  working  out  the  proper  relationships.  Since  fhe  nace''es  are  not 
generally  limiting  components,  a  side  bending  allowable  <s  considered 
reasonable. 
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DIFcRArTinN  PHASF  LOADF 


4.0 

The  question  investigated  by  this  study  is  "What  is  the  importance  of  inr'udmq 
the  diffraction  phase  loads  in  calculations  of  structural  response  fo  a  nw  'oar 
gust'’"  The  VIBRA-4  program  was  chosen  for  the  comparison  work.  It  was  ’he  or  y 
program  investigated  which  met  the  needs  of  the  study.  The  NOVA-! 1  orogram 
includes  the  diffraction  phase  loads  but  not  the  "long"  time  response  of  maior 
structural  components.  The  VIRRA-6  program  currently  does  not  have  the 
diffraction  phase  loads  incorporated  although  the  equations  for  them  are 
included  in  Appendix  A  of  Reference  <3. 

The  VIBRA-4  program  was  easily  changed  to  bypass  the  diffraction  phase  loads, 
the  aerodynamic  loads  described  in  Section  i.7  wer«  modified  in  thp  following 
way.  Thp  part  of  the  program  that  tests  for  the  non-(t->mensional  time 
S-intnersion  and  calculates  the  diffractive  phase  loads  was  bypassed  and  thp 
resulting  normal  force  coefficient  then  used  is  as  shown  in  Figure  S.  This  can 
be  compared  to  Figure  ?  which  shows  the  normal  force  coefficient  use-i  with  the 
diffraction  phase  loading  included. 

Two  different  airplanes  were  analyzed  using  the  version  of  VIBRA-4  described  in 
Section  3.0.  The  mathematical  model  of  thp  KC-135  used  for  the  study  is 
described  in  Section  4.1.  The  mathematical  model  of  the  747-200  used  as  the 
other  aircraft  in  the  study  is  described  in  Section  4 .2. 

Three  blast  orientations  were  used  for  the  study;  vert ical -bpneat h , 
vertical-overhead,  and  lateral-from  the  right  side  of  the  vehicle.  the  same 
input  data  decks  were  used  for  the  analysis  without  the  diffraction  phase  7oads 
as  described  above. 

The  comparisons  between  the  VIBRA-4  results  with  diffraction  phasp  loa^s  and  the 
VIBRA-4  results  without  diffraction  phase  loads  are  shown  ’n  sect’on  4. a. 

4.1  Oynamic  Model  -  KC-P5 

The  mathematical  model  for  the  Kf,-t75  used  to  conduct  the  present  analysis  was 
previously  developed  and  documented  in  Reference  ?.  The  model  consists  of  a 
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discrete  mass  representation  o*  the  aircraft  inertia  properties  by  l'7  total 
masses;  28  masses  on  each  wing,  8  of  which  represent  the  two  engines,  7ri 

fuselage  masses,  14  masses  on  each  horizontal  stabilizer  and  19  masses  on  the 
vertical  stabilizer.  The  total  gross  weight  represented  by  these  masses  ’s 

292,000  lbs.  The  distribution  of  masses  was  selected  to  represent  each  elastic 
axis  panel  weight,  center  of  gravity  and  torsional  inertia  about  the  elastic 
axis.  The  only  exceptions  were  the  engine  nacelles  for  which  both  pitch  and 
roll  inertia  were  matched.  The  wing  masses  were  initially  in  the  wing  chord 
plane  since  the  VIBRA-4  program  requires  mass  points  associated  with  the 
aerodynamic  panels  to  be  at  the  same  z  coordinate.  This  simplification  reduces 
the  model  complexity  without  significantly  affecting  the  vihr?t-jon 
characteristics.  Tn  addition,  the  model  consists  of  15  syimietric  and  15 

antisymmetric  free  aircraft  coupled  modes  of  vibration.  The  coefficients  of 
these  modes  are  used  in  Equation  19  to  determine  the  response  of  the  structure. 
These  modes  were  derived  from  vibrational  analysis  based  upon  properties  formed 
from  the  discrete  mass  representation  and  the  stiffness  properties  of  the 
structure. 

In  addition  to  the  inertia  and  modal  data,  the  mathematical  model  includes 
aerodynamic  data.  In  these  data,  the  lifting  surfaces  and  the  fuselage  are 

divided  into  aerodynamic  strips  or  panels.  There  is  a  total  of  45  aerodynamic 
panels  on  the  complete  aircraft  with  11  on  each  wing,  4  on  the  fuselage,  6  on 
each  horizontal  stabilizer  and  7  on  the  vertical  stabilizer.  An  area  and  its 
centroidal  coordinates  are  associated  with  each  panel. 

Additional  components  of  the  mathematical  model  include  bending  moment  da*a  for 
each  component,  and  engine  thrust  data.  The  bending  moment  data  consists  of  *-he 
coordinates  of  the  point  on  the  elastic  axis  of  the  component  where  bendinq 
calculations  are  required,  and  the  limiting  bending  moment  at  that  station. 
This  limiting  bending  moment  can  either  be  the  ultimate  or  catastrophic  damage 
load  or  it  can  be  the  safe,  or  no  damage  load.  In  this  way,  the  response  for 
any  desired  loading  allowable  can  be  obtained.  To  assist  m  'ocatinq  stations 
on  each  component.  Table  1  provides  a  cross  reference  between  buttock  line, 
station  and  percent  semispan  for  each  of  the  bending  moment  stations.  The 
actual  loads  used  as  allowables  are  shown  in  Table  2. 
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TABLF  ' 

KC-135  COMPONENT  BENDING  MOMENT  LOCATIONS  FOR  GIJST  ANALYSIS 


COMPONENT 

STA. 

NO. 

STATION 

n 

Wing 

1 

WBL  157. ’ 

.20 

? 

WBL  259.6 

.33 

3 

WBL  t??.5 

.409 

4 

WBL  S03. 5 

.64 

5 

WBL  650.7 

.^01 

Stabil i zer 

1 

SBL  41.2 

.17’ 

2 

SBL  BO.O 

,T?3 

0 

SBL  119.0 

.dQK 

Fin 

1 

5WL  353.3 

.204 

.  2 

RWL  402.5 

.  34  P 

Body 

1 

BS  inRO 

Nacelles 

2 

3 

1 

BS  1200 

BS  1360 

BWl  1  74. 714 

WRL  -  Wing  Ruttock  Line 

BS  - 

Body  station 

SBL  -  Stabilizer 

Buttock 

Line  BWL 

-  Body  Waterl  i 

H  -  Percent  Semispan 


r  /-  '  -  r  ~  - ~rr 


TABLE  2 

KC-135A  ALLOWABLE  LOADS  FOR  MAJOR  COMPONENTS 


COMPONENT 
W1  ng 


Fusel  age 


Hori zontal 
Stabilizer 


Vertical 
Stabll Izer 


Engine 

Nacelle 


ULTIMATE 

1001  LIMIT 

STA. 

LOAD 

.LOAD 

NO. 

STATION 

106  In 

-lb 

10°  In-lb 

1 

WBL  157.3 

82.7 

Up 

51.13 

53.6 

Dn 

35.73 

2 

UBL  259.6 

51.4 

Up 

34.27 

34.5 

Dn 

23.00 

3 

WBL  503.5 

15.8 

Up 

26.60 

27.3 

Dn 

18.20 

4 

WBL  503.5 

15.8 

Up 

10.53 

10.3 

Dn 

6.87 

5 

WBL  550.7 

12.0 

Up 

8.0 

7.8 

Dn 

5.20 

1 

BS  1080 

43.3 

Up 

28.87 

74.6 

Dn 

49.73 

52.6 

Side 

35.07 

2 

BS  1200 

32.1 

Up 

21.4 

47.3 

Dn 

31.53 

44.7 

Side 

29.8 

3 

BS  1360 

19.8 

Up 

13.2 

27.4 

Dn 

18.27 

29.1 

Side 

19.4 

1 

SBL  41.23 

3.47 

l*> 

2.31 

3.864 

Dn 

2.576 

2 

S8L  80.0 

1.875 

Up 

1.25 

2.333 

Dn 

1.555 

3 

S8L  119.0 

1.12 

Up 

0.747 

1.414 

Dn 

0.943 

1 

BWL  353.3 

4.11 

Side 

2.74 

2 

BWL  402.5 

2.34 

Side 

1.56 

Inbd. 

BWL  174.7 

0.73 

Side 

0.487 

(Xitbd. 

BWL  174.7 

0.666 

Side 

0.444 

WBL  -  Wing  Buttock  Line  SBL  -  Stabilizer  Buttock  Line 

BS  -  Body  Station  WL  -  Water  Line 
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4.2  DYNAMIC  MODEL  -  747-200 


The  mathematical  model  for  the  747-200  used  to  conduct  the  present  analysis  was 
previously  developed  and  docunented  in  Reference  7.  The  mode1  consists  of  a 
discrete  mass  representation  of  the  aircraft  inertia  properties  by  177  tota1 
masses;  49  masses  on  each  wing,  0  of  which  represent  the  two  engines,  74 
fuselage  masses,  15  masses  on  each  horizontal  stabilizer  and  15  masses  on  the 
vertical  stabilizer.  The  total  gross  weight  represented  by  these  masses  is 
7Q5,000  lbs.  The  distribution  of  masses  was  selected  to  represent  each  elastic 
axis  pane?  weight,  center  of  gravity  and  torsional  inertia  about  the  elastic 
axis.  The  only  exceptions  were  the  engine  nacelles  for  which  both  pitch  and 
roll  inertia  were  matched.  The  wing  masses  were  initially  in  the  wing  chord 
plane  since  the  VIBRA-4  program  requires  mass  points  associated  with  the 

aerodynamic  panels  to  be  at  the  same  z  coordinate.  This  simplification  reduces 
the  model  complexity  without  significantly  affecting  +  he  vibration 
characteristics.  In  addition,  the  model  consists  of  1?  symmetric  and  1? 
antisymmetric  free  aircraft  coupled  modes  of  vibration.  The  coefficients  of 
these  modes  are  used  in  Equation  19  to  determine  the  response  of  ^he  structure. 
These  modes  were  derived  from  vibrational  analysis  based  upon  properties  formed 
from  the  discrete  mass  representation  and  the  stiffness  properties  of  +he 
structure. 

In  addition  to  the  inertia  and  modal  data,  the  mathematical  model  includes 

aerodynamic  data.  In  these  data,  the  lifting  surfaces  and  the  fuselage  are 

divided  into  aerodynamic  strips  or  panels.  There  is  a  total  of  4F>  aerodynamic 
panels  on  the  complete  aircraft  with  17  on  each  wing,  4  on  the  fuselage,  5  on 
each  horizontal  stabilizer  and  5  on  the  vertical  stabilizer.  An  area  and  its 
centroidal  coordinates  are  associated  with  each  panel. 

Additional  components  of  the  mathematical  model  include  bending  moment  data  for 
each  component,  and  engine  thrust  data.  The  bending  moment  data  consists  of  the 
coordinates  of  the  point  on  the  elastic  axis  of  the  component  where  bending 

calculations  are  required,  and  the  limiting  bending  moment  at  that  station. 
This  limiting  bending  moment  can  either  be  the  ultimate  (or  catastrophic!  damage 
load  or  it  can  be  the  safe,  or  no  damage  load.  In  this  way,  the  response  for 
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TABLE  3 

747-200  COMPONENT  BENDING  MOMENT  LOCATIONS  FOR  GUST  ANALYSIS 


STA. 


COMPONENT 

Wing 

NO. 

1 

STATION 

WBL  176.1 

n 

.150 

2 

WBL  419.1 

.357 

3 

WBL  584.7 

.498 

4 

WBL  743.1 

.633 

5 

WBL  880.5 

.750 

Stab i 1  i  zer 

1 

SBL  100.0 

.229 

2 

SBL  231.8 

.530 

Fin 

1 

BWL  385.15 

.052 

2 

BWL  537.05 

.443 

Body 

1 

BS  1640. 

2 

BS  2180. 

Nacelles 

1 

BWL  199.9 

-  Wing  Buttock  Line 

-  Body  Station 

-  Percent  Semi  span 


SBL  -  Stabilizer  Buttock  Line 
BWL  -  Body  Water  Line 
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TABLE  4 


747-200  ALLOW  ARL:  LOADS  FOR  MAJOR  COMPONENTS 


ULTIMATE 

100*  LIMIT 

STA. 

LOAD 

10°  in-lb 

LOAD 

COMPONENT 

NO. 

STATION 

ID6  in-lb 

Wing 

1 

WRL 

1  76.1 

109. 5 

Up 

206.il 

281.86 

Dn 

189.24 

2 

WBL 

41°. 1 

1 79 , 1 5 

Up 

1 1 9.41 

141.56 

nn 

14 . 71 

7 

WP' 

6P4.  ’ 

108. 1° 

"p 

>h 

27.26 

Dn 

IP. .  ' 

4 

WB1 

M7.  1 

56.84 

Up 

i-’.pg 

26.14 

Dn 

>4. no 

5 

WRl 

p«o.  ^ 

24.19 

Up 

16.11 

15.16 

Dn 

10.7* 

F uspI age 

1 640. 

115.0 

Up 

on.  ap 

120.0 

Dn 

7\7  77 

49.n 

c’de 

1 2 1  6  7 

Hori zonta1 

1 

SBL 

mo.n 

14.o 

Up 

0.91 

Stabii i zer 

19.1 

on 

1?.  77 

2 

SBL 

?il.R 

4.6 

Up 

7.P7 

6.0 

Dn 

4.  op 

Vertu  al 

Stab ’I  izer 

1 

ML 

185.15 

11.6 

Side 

I.O7 

2 

WL 

5.17.05 

7.4? 

Side 

4.95 

Engine 

Nac  pi  1  p 

Tnbd. 

WL 

199.0 

90.2 

Side 

60.2 

Outbd. 

WL 

199.7 

12.5 

Side 

48.li 

WBL 

-  Wing  Buttock  Line 

SBL 

-  Stabilizer  Buttock  Line 

BS 

-  Body  Station 

WL 

-  Water  Line 

any  desired  loading  a’lowable  can  be  obtained.  To  assist  ’orating  stat'ons 
on  each  component.  Table  3  provides  a  cross  reference  between  buttock.  'ine, 
station  and  percent  semispan  for  each  of  the  bending  moment  stations.  The 
actual  loads  used  as  allowables  are  shown  in  Table  4. 

4.3  COMPARISON  OF  VIBRA-4  RESULTS 

The  VIBRA-4  output  for  each  run  included  the  bending  moment  at  several  stations 
for  each  time  increment.  At  the  end  of  each  time  history  the  ratio  of  the  peat 
load  to  the  allowable  load  at  each,  station  is  also  printed  out. 

Table  5  shows  the  flight  condition  used  tor  the  KC-P5  analysis.  The  time 
histories  are  shown  in  Figures  A  through  1°.  Most  of  these  show  small 
differences  between  the  VIBRA-4  results  and  the  VIRRA-4  without  diffraction 
phase  loads.  Table  6  is  a  summary  of  the  peak  bending  moment  ratios  calculated 
for  the  KC-P5.  The  summary  shows  the  VIRPA-4  ratios,  without  diffraction  phase 
loads  ratios,  and  the  percent  increase  in  the  ratio  by  including  the  diffraction 
phase  loading. 

TABLE  5 

KC-135  FLIGHT  CONDITION 


Vel ocity 
Altitude 
Overpressure 
Yield 

Material  Velocity 
Range 


729  fps-true 
20,000  feet 
1.07  psi 
1.0  MT 

no. 06  fps 

25,000  fee* 


The  747-200  flight  condition  analyzed  is  shown  in  Table  7.  The  t>me  histories 
from  this  analysis  are  shown  in  Figures  20  through  38.  Like  the  KC-P5  analysis 
some  show  only  small  differences  between  the  VIBRA-4  results  and  the  VTBRA-4 
without  diffraction  phase  loads.  However  the  time  histories  for  other  stations 
show  vast  differences  in  peak  loads  between  the  two  sets  of  results.  The 
sunmary  of  peak  bending  moment  ratios  in  Table  8  show  that  increases  in  peak 
loads  of  300  percent  were  calculated  as  a  result  of  including  the  diffraction 
phase  loads. 
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TARLF.  6 


KC - 1 35  PEAK  BENDING  MOMENT  RAT  I n  SUMMARY 


STA.  BLAST  FROM  BELOW 
n  1  ?  t 


Right  Wing 

.20 

.33 

.409 

.64 

.701 

Left  Wing 

.20 
.33 
.409 
.64 
.  701 

Fusel  age 

1080 

Station 

1200 

1360 

Right  .17 

Horizontal  .SB 

Stabilizer  .50 


Left  .17 

Horizontal  .33 

Stabilizer  .50 


Fin  .204 

.348 


539 

.520 

+3.6 

616 

.598 

+  3.0 

585 

.678 

+  1.2 

627 

.582 

+  7.7 

563 

.509 

+  10.6 

.494  .  '35R  +38.0 

.627  .450  +39.3 

.630  .448  +40.6 


.360  .311  +15.7 

.435  .372  +16.9 

.473  .376  +25.8 


BLAST  FROM  ABOVE 
1  2  3 

. 2°1  .261  +3.8 

.302  .298  1.3 

.291  .289  0.7 

.280  .253  10.7 

.  24'!  .226  3.5 


.813 

.645 

26.0 

.886 

.687 

28.9 

.841 

.636 

32.2 

378 

.398 

-5.0 

574 

.516 

+11.2 

733 

.565 

+29.  v 

BLAST 

FROM  RIGHT 

1 

2 

.276 

.269 

+2.6 

tj-j 

#  ■ 

.307 

8.6 

.323 

.302 

6.9 

.302 

.239 

8.2 

.276 

.256 

7.8 

.293 

.287 

2.1 

.340 

.325 

3.4 

.336 

.313 

7.3 

.282 

.262 

7.6 

.266 

.238 

11.8 

.352 

.269 

35.9 

.367 

.275 

33.5 

.335 

.254 

31.9 

.340 

.271 

25.5 

.348 

.  2°6 

17.6 

.375 

.281 

19.2 

.297 

.271 

9.6 

.327 

.'’86 

14.3 

.315 

.265 

18.9 

.3?8 

.577 

28.8 

.922 

.699 

71.9 

COLUMN  1  VIBRA-4 

COLUMN  2  WITHOUT  DIFFRACTION  PHASE  LOADS 

COLUMN  3  [COLUMN  1/COLUMN  2  -  1.03  X  100% 
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Because  the  user  cannot  know  in  advance  the  effect  the  diffraction  phase  loads 
will  have  on  the  total  structural  response  of  his  system  or  on  the  magnitude  of 
change  in  peak  load  at  any  one  station,  the  diffraction  phase  loads  should  be 
included  in  the  structural  response  calculations  for  the  system. 


747-2(10 

Vel ocity 
A1 1  i  t  ude 
Overpressure 
Yiel  d 
Material  Velocity 
Range 


TABLE  7 

FLIGHT  CONDITION 

674  fps-true 
25000  feet 
1  .*7  psi 
5.0  MT 
107.9  fps 
70000  feet 
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table  p 


747-200  PEAK  RENDING  MOMENT  RATIO  SUMMARY 


STA. 

BLAST 

FROM 

BELOW 

RLAST 

FROM 

ABOVE 

BLAST 

FROM 

RTGHT 

n 

2 

1 

2 

1 

2 

■> 

Right  Wing 

.150 

.482 

.417 

15.2 

.707 

.577 

-47.1 

.742 

.728 

4.7 

.757 

.456 

.408 

11.0 

.717 

.660 

-47.4 

.747 

.727 

4.7 

.TV 

.4/7 

,  447 

7.7 

.260 

.46Q 

-44.4 

.  7R2 

.776 

17.8 

.  750 

.pig 

.460 

78.1 

.  ?8Q 

.707 

-27.7 

.  4no 

.747 

10.7 

Loft  Wing 

.  '  S° 

it  ’ 

.40'- 

k'- 

.421 
.428 
,  746 
.IMS 

.406 

.707 

.774 

f117 

4.0 

7.0 

7.2 

14.4 

Fusel  age 

St  at  ion 

1  640 

2180 

.28’ 

1  .017 

.761  -20.6 
.452  +124.0 

.572 

1  .187 

.450 

.775 

+27.1 

+57.1 

1  .  783 
1.847 

1  .077 

.999 

66.2 

84. 

Right 

Hori zontal 

St  ah i 7 i zer 

.229 

.530 

.987 

1.204 

.661  57.1 

.790  52.7 

1  .200 
1.768 

.774 

.853 

54.9 

60.2 

.616 

.679 

.244 

.242 

152.1 

180.1 

Left 

Hori zontal 
Stabil izer 

.229 

.530 

.804 

.8.94 

.227  250.1 
.219  280.2 

Fi  n 

.052 

.443 

1.842 

2.6QR 

.442 

.641 

336.1 

704.8 

COLUMN  1 

VIBRA-4 

COLUMN  2 

WITHOUT  DIFFRACTION  PHASE 

LOADS 

COLUMN  3 

[COLUMN  1 /COLUMN  2  -  1.0] 

x  100% 
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5.0 


CONCLUSIONS 


The  results  in  Section  4.0  show  that  the  diffraction  phase  loads  should  bn 
included  in  '■.he  structural  response  calculations.  The  magnitude  of  the 
diffraction  phase  1  .-ads  was  not  large  enough  to  change  the  peak  loads  alone,  but 
the  diffraction  phase  loads  added  to  the  gust  loads  can  have  a  significant 
effect  on  the  peak  loads  calculated.  This  effect  varies  from  station  to  station 
and  from  one  aircraft  to  another.  There  is  no  known  rule  of  thumb  or  equations 
that  would  help  the  user  determine  the  magnitude  of  the  effect  of  including 
diffraction  phase  loads  in  the  nuclear  gust  response  calculations.  Therefore 
structural  response  codes  such  as  VIBRA-4  and  VIRRA-6  should  include  the 
diffraction  phase  loading  with  the  gust  response  loading. 


6.0 


RECOWENDATIQNS 


1.  VIBRA-4  should  not  be  modified  to  remove  the  diffraction  phase  loads. 

2.  VIBRA-6  should  have  the  diffraction  phase  loads  added. 

3.  Further  studies  should  be  made  to  determine  the  effect  the  diffraction 
phase  loads  have  on  the  peak  loads  with  variations  in  altitude  for  a 
given  aircraft  at  various  speeds  and  vice  versa. 

4.  Additional  studies  should  be  made  to  determine  the  effect  the  diffraction 
phase  loads  have  on  the  peak  loads  with  size  of  the  vehicle.  Small 

air  vehicles  such  as  a  swept  wing  fighter  and  a  winged  missile  such  as 
the  Air  Launched  Cruise  Missile  (ALCM)  could  be  used  for  the  input  models. 
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Figure  I  Gust  analysis  coordinate  system. 
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function  of  aerodynamic  time 


VIEW  LOOKING  AFT 


+  X 

NOTE  -  AXIS  SYSTEM  IS  V I  BRA-4  LOCAL  AXIS  SYSTEM 


Figure  3  Fuselage  allowable  moment  boundary. 
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Figuri  4  Overpressure  as  a  function  of  range  for  a  l-KT  explosion  at  sea  level. 
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Figure  7  KC-T35  right  hand  wing  -  blast  from  below 
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KC-135  body  vertical  bending  -  blast  from  below 
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Figure  12  KC-135  right  hand  horizontal  stabilizer  -  blast  from  above 
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Figure  16  KC-135  left  hand  horizontal  stabilizer  -  blast  *rom  the  side 
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Figure  17  KC-135  vertical  bocr/  benc'rg  -  blast  from  the  side 


Figure  18  KC-135  horizontal  oca-'  a?ng  -  blast  from  the  side 
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•ENDING  MOMENT 


Fig or*  21  747  rfghf  hand  wing  -  blasf  from  oeiow 
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Figure  22  747  right  hand  wing  -  blast  from  below 
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Figure  25  747  vertical  body  bending  -  blast  from  below 
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Figure  27  747  right  hand  wing  -  blast  from  above 


Figure  28  747  right  hand  wing  -  blast  from  above 
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Figure 


•ENDING  MOMENT 
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Figure  30  747  right  hand  horizontal  stabilizer  -  blast  from  above 
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Rgur«  31  747  vertical  body  bending  -  blast  from  above 
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Figure  32  747  vertical  body  bending  -  blast  from  above 
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Figure  33  747  right  hand  wing  -  blast  from  side 
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Figure  34  747  right  hand  wing  -  blast  from  side 
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Figure  35  747  left  hand  horizontal  stabilizer  -  olasr  rrem  »iae 
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Figure  37  747  horizontal  body  bending  -  blast  from  side 
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Figure  38  747  horizontal  body  bending  -  blast  from  side 
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Figur*  39  747  fin  -  blast  from  side 
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